Attenuation of shear sound waves in jammed solids 
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We study the attenuation of long-wavelength shear sound waves propagating through model 
jammed packings of frictionless soft spheres interacting with repulsive springs. The elastic attenua- 
tion coefficient, a(u>), of transverse phonons of low frequency, u, exhibits power law scaling as the 
packing fraction <j> is lowered towards <j> c , the critical packing fraction below which rigidity is lost. The 
elastic attenuation coefficient is inversely proportional to the scattering mean free path and follows 
Rayleigh law with a(uj) ~ oj 4 ((f> — <^>c) _5//2 for u much less than oj* ~ (<j> — C ) 1 ' /2 , the characteristic 
frequency scale above which the energy diffusivity and density of states plateau. This scaling of the 
attenuation coefficient, consistent with numerics, is obtained by assuming that a jammed packing 
can be viewed as a mosaic composed of domains whose characteristic size £* ~ (<f> — <^ c ) _1/ ' 2 diverges 
at the transition. 

PACS numbers: 45.70.-n, 61.43.Fs, 65.60.+a, 83.80.Fg 



I. INTRODUCTION 

The challenge of attaining a robust understanding of 
the elasticity and vibrational dynamics of amorphous 
solids continues to be an elusive goal in material science 
3]. One of the main difficulties arises from the uncer- 
tainty surrounding the application of the time honored 
formalism of continuum mechanics to disordered struc- 
tures characterized by strong spatial inhomogeneities and 
by a non affine response to external perturbations, which 
often acts as a precursor to elastic instabilities (j-Q . The 
limitations of our current understanding are manifest 
when one considers the slow pace of progress in achieving 
a satisfying picture of the low temperature properties of 
glasses (such as the ubiquitous presence of a plateau in 
their thermal conductivity) or of sound propagation in 
disordered solids such as emulsions, foams and granular 
media 041 EMU ■ 

Several authors have been intrigued by the prospect of 
a unified description of these apparently diverse phenom- 
ena that relies on general principles common to molecu- 
lar glasses as well as amorphous packings of bubbles or 
grains 0, [2(| H3|- The generic starting point, often in- 
voked in the context of simple models and computer sim- 
ulations, is the existence of a characteristic mesoscopic 
length scale, intermediate between the microscopic size 
of the basic components and the sample size, below which 
continuum elastic theories cease to be valid [1, [HI ■ 
A consequence of the breakdown of classical elasticity is 
that the vibrational degrees of freedom characteristic of 
amorphous solids cannot always be identified as simple 
plane waves. Thus, mechanical and transport proper- 
ties controlled by the vibrational dynamics of amorphous 
solids differ markedly from their crystalline counterparts. 

Experimental progress in measuring elastic wave scat- 
tering has provided reliable nondestructive methods to 
probe the structure of amorphous materials at different 
length scales [22| . By varying the frequency or the wave- 
length of the incident wave, different spectral regimes of 
the vibrational dynamics of the disordered samples be- 
come experimentally accessible. This provides an iden- 



tikit of the vibrational degrees of freedom of the amor- 
phous solids under investigation and paves the way to un- 
derstanding how they emerge from characteristic length 
scales inherent to these disordered structures. 

In this work we concentrate on how sound propagates 
through jammed packings of frictionless soft repulsive 
spheres whose density can be varied above the critical 
threshold necessary to insure mechanical rigidity [24j . 
This simple and concrete model system can be readily 
generated on a computer and captures several of the 
generic properties of amorphous materials, such as the 
onset of excess vibrational modes at a characteristic fre- 
quency uj* [25j . above which the scattering mean free 
path and the wavelength are comparable [l|, The 
characteristic frequency lo* heralds a spectral regime of 
strong scattering characterized by a small and frequency 
independent energy diffusivity. This feature is also seen 
in experiments carried out on glasses as well as porous 
media [HHl]. 

A unique and intriguing feature that motivates our 
choice of jammed solids over other model systems is the 
fact that their elastic properties exhibit scaling behavior 
with packing fraction or connectivity above the critical 
density <fi = <f> c (Point J) at which the particles are just 
touching each other [24| . In particular, the mesoscopic 
length I* ~ (<j> — <\> c )~ x ' 2 above which continuum elas- 
ticity is expected to hold in this model, can be made 
to diverge as the unjamming transition is approached by 
progres sively decompressing the sample towards point J 
20, 21]. However, some uncertainty still exists regarding 
the existence and nature of this diverging length scale. 

A recent effective-medium investigation of phononic 
transport in random networks of harmonic springs cap- 
tured several properties shared by jammed packings of 
soft spheres [32|]. The random network model predicts 
a sharp crossover at w* in the energy diffusivity from a 
plateau to a low-o; Rayleigh scattering regime, with a 
pre-factor that scales with connectivity and eventually 
vanishes at the isostatic point. Despite its success, this 
mean field approach does not yield any transport signa- 
tures of the characteristic length scale £* that controls 
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the fluctuations in the microscopic structure. 

The diffusivity plateau above ui* was also studied di- 
rectly for the jammed packings by evaluating the Kubo 
formula without recourse to a mean field approximation 
[l], 0. However, no attempt was made in Ref. 0, EI to 
extract the packing fraction dependence of Rayleigh law 
for lo -C uj*. 

In this article, the scaling of the pre-factor to Rayleigh 
law is determined by studying the attenuation coefficient 
of long-wavelength transverse acoustic phonons which 
carries a vestige of the important length scale £* asso- 
ciated with the unjamming transition. The numerical 
data of Ref. [l|, E| is re-analyzed and shown to be consis- 
tent with the present analysis, but not with mean field 
calculations because the latter do not yield the length 
scale £*. Our scaling analysis of the attenuation coeffi- 
cient provides insights on the structure of the jammed 
solids. In particular, it suggests a simple view of the 
jammed solid as a mosaic composed of domains whose 
characteristic size is of order £* . 

In addition to its relevance as a probe of the mi- 
crostructure of a jammed solid, the present analysis sug- 
gests a methodology to infer the fabric of a granular me- 
dia from experimental investigations of sound propaga- 
tion under external loading 

II. REVIEW OF THE ELASTICITY OF 
JAMMED SPHERE PACKINGS 

The model jammed solids studied in this work are 
amorphous packings of frictionless spheres generated by 
conjugate-gradient energy minimization using the proto- 
col described in Ref. |24j |. The harmonic pair potential 
V(rij) between particles is purely repulsive: 

V ( r ij) = ~(X ~ nj / (Tijf if r lJ <a ij 

V(rij) = if r v] > , (1) 

where r$j denotes the distance between the centers of 
particles i and j and <xy is the sum of their radii. 

At the unjamming transition (j) — (j) ci the average coor- 
dination number, z, increases discontinuously from zero 
to the "isostatic" value z c = 2D, where D is the dimen- 
sionality of the sample 29]. Above <p c at A(f> = 4> — (f> c , 
the average coordination number exhibits power law scal- 
ing [Mill: 

Az = z - z c ~ Ac/) 1/2 . (2) 
In addition, the shear modulus scales with A<p according 

to mm 

G (3) 

whereas the bulk modulus B is only weakly dependent on 
packing fraction near the jamming point. The anomalous 
scaling of the elastic moduli reflects the importance of 
non-affine deformations |3l| . 



The transverse and longitudinal speeds of sound Vt and 
vi are proportional to the square root of the shear and 
bulk moduli, which implies that only the transverse speed 
of sound scales with A(j> according to 

1H (4) 

Upon substituting Eq. @ into the Deybe formula, the 
ratio of the transverse to the longitudinal density of states 
at low frequency diverges as A<p — >• 0. Above a character- 

1 /2 

istic frequency uj* ~ A(f> ' the density of states exhibits 
a plateau that extends all the way down to zero frequency 
at point J [25|. That is to say, in a jammed solid right 
at the threshold of rigidity the familiar Deybe regime 
is completely absent and the acoustic waves, which are 
the familiar excitations of continuum elastic media, are 
suppressed. 

An explanation of this puzzling phenomena was pro- 
posed by Wyart and collaborators who have suggested 
that a characteristic length scale £*, which diverges as 
4> appr oaches cf> c , controls the elasticity of jammed solids 
[2Cll2f| . Here, we briefly review their "cutting argument" 
originally proposed to obtain the packing fraction depen- 
dence of £* . 

Consider the mechanical stability of a jammed packing 
at a packing fraction <fi > <f> c , such that there is a den- 
sity of inter-particle contacts Az in excess of the mini- 
mum isostatic value z c necessary to ensure stability. The 
properties of a continuum elastic medium should remain 
invariant upon cutting it in two parts. If you cut a spher- 
ical blob of radius £ out of a jammed medium, a number 
of contacts equal to the area £ d ~ 1 of the blob will be bro- 
ken. However, the excess number of contacts that exist 
in the blob is given by £ d Az. One obtains the size £* of 
the smallest blob that will remain rigid upon cutting, by 
balancing the number of broken contacts with the excess 
ones. This leads to [H EH 



Az ' w 

where Eq. ([2]) was used to obtain the explicit dependence 
on packing fraction. According to this elegant argument, 
below the length scale £* the medium is to be considered 
floppy and hence unable to support the acoustic (quasi) 
plane waves that give rise to the Deybe counting in the 
vibrational density of states. Long wavelength phonons 
with A 3> £* can still be supported because on those 
length scales the jammed packing is rigid and can be 
described as a continuum elastic medium. 

The conclusions of this simple argument are consistent 
with the numerical studies carried out by Ellenbroek and 
collaborators who probed the fluctuations in the force re- 
sponse of a jammed solid to inflating the diameter of a 
single particle at the center of the sample |3l| . The spa- 
tial range of these fluctuations is given by £* and grows 
upon approaching the jamming point. 

Despite this progress, there is still uncertainty regard- 
ing the exact character and the very existence of this im- 
portant length scale. In particular it is not clear whether 
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£* can be viewed as a structural signature of the un- 
jamming transition or merely as a dynamic length scale 
that appears in the vibrational dynamics of these fragile 
solids. 

In the next section, we will address this issue by study- 
ing the packing fraction dependence of the scattering at- 
tenuation coefficient of long wavelength phonons. Such 
waves can be excited by introducing in the jammed sys- 
tem an oscillatory perturbation of very low frequency to 
probe dynamically the divergent fluctuations observed in 
Ref. (Hi. 



III. PACKING FRACTION DEPENDENCE OF 
THE ATTENUATION COEFFICIENT 

A. Scaling analysis. 

Consider a shear sound wave of very low frequency 
lu propagating through a jammed packing with a well 
defined transverse speed of sound, v t . In this section, we 
develop simple physical insights that allow us to infer the 
presence of the diverging length scale £* from the elastic 
scattering of these long wavelength shear waves. 

The approach pursued here is different from previous 
studies which have sought signatures of this length scale 
in the vibrational dynamics around and above lu*. The 
regime lu ^ lu* has generally been considered uninter- 
esting because the density of states of weakly scattered 
phonons n(w) is obtained by simply counting how many 
waves you can fit in a three dimensional box. As point 
J is approached, the shear modulus becomes much softer 
than the bulk modulus (see Eq. §5§ and the follow- 
ing comment) and the spectrum is dominated by shear 
waves. The transverse density of states nt(w) takes up 
the generic form 



n t (w) 



(G) 



The packing fraction dependence of the pre-factor of 
Deybe law is entirely determined by the scaling of the 
shear modulus in Eq. ([3]) and reads 



n t (w) 



A</>' 



3/4 



(7) 



Equation ([6J is merely what you expect for a perfectly 
ordered solid, hence it does not carry distinctive signa- 
tures of the amorphous structure of a jammed packing. 
This has prompted various investigators to seek indirect 
evidence for £* at higher lu where the density of states 
departs dramatically from Deybe behavior and exhibits 
a characteristic plateau. 

Here we show that the sound attenuation coefficient 
a(uj) is a more sensitive spectral measure than the den- 
sity of states and can display anomalies as a result of 
elastic scattering from structural inhomogeneities even 



in the low lu regime for which the density of states ap- 
proximately follows Deybe law. 

In order to define the sound attenuation coefficient, 
denote by Iq the initial intensity (proportional to the 
elastic energy density) and by I(x) the reduced intensity 
attained by the wave after traveling a distance x through 
the sample. The coefficient of sound attenuation a{ui) 
can then be extracted from the relation 



I(x) = I e- aiu])x 



(8) 



In what follows we focus our attention to the simpler 
case in which sound attenuation is primarily triggered 
by elastic scattering. This restriction precludes us from 
considering the effects of absorption and dissipation due 
to friction which can be significant for a real granular 
packing [Hp. 

The elastic scattering contribution to the attenua- 
tion coefficient is inversely proportional to the scattering 
mean free path £ s (lu) (26j 



a(w) 



1 



(9) 



Sound attenuation in granular media was observed to be 
dominated by transverse waves [ljj, [ll|. According to 
Rayleigh's law, the mean free path £ s (lu) of a tranverse 
sound wave is given by the generic form [H, H3] 



■<»>-(=)' 



(10) 



where D is the length scale that characterizes the disor- 
dered structures responsible for the scattering. As an- 
ticipated, the presence in Eq. (fTUf of the multiplicative 
factor D~ 3 , which does not appear in Eq. ([5]), makes 
the mean free path £ s (lu) a more sensitive probe of amor- 
phous structures than the density of states n t {uu). 

In polycrystalline solids, the length D is given by the 
typical distance between grain boundaries. If the char- 
acteristic length scale of structural inhomogeneities in a 
jammed media had no packing fraction dependence, Eq. 



(|TU| would lead to t 3 (u) 



A(f>. Instead, we make 



the working assumption that the divergent fluctuations 
of size £* ~ (4> — 'Ac)" 1 / 2 : observed in response to a static 
perturbation in Ref. [3l|, are still crucial when the sys- 
tem is probed dynamically by sending low lu elastic waves 
through it. 

Upon setting the speed of sound equal to the square 
root of the shear modulus and D ~ £* in Eq. (fTO]) , we 
obtain 



1 G 2 



(11) 



Substituting Equations (|3]) and ([5]) into Eq. (TTTj) gives 
the packing fraction dependence of the attenuation coef- 
ficient (or equivalently the inverse of the mean free path) 



a(LU) 



(12) 
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Equation (|T2|) suggests that in the weak scattering regime 
u) <C U3* , the absorption length l/a(uj) vanishes as point 
J is approached in a way controlled by the anomalous 
scaling of the shear modulus and of the length scale I* , 
providing indirect evidence for the existence of the latter. 
Indeed, a simple physical picture of how elastic waves 
propagate through a jammed packing under loading is 
obtained from viewing this disordered solid as a mosaic 
of domains with characteristic size I* ~ (4> — </> c )~ 1//2 . 
The resulting mean free path £ s is not simply given by the 
characteristic length I* . Instead, upon approaching point 
J, £ s decreases as I* ~ A0" 1 / 2 increases, as described by 
Eq. (1111) . This is consistent with the intuitive notion that 
the unjamming transition is accompanied by a gradual 
loss of rigidity that decreases the efficiency of the energy 
transport by plane wave excitations hence resulting in a 
vanishing pre- factor for Rayleigh law at point J. 

The noteworthy feature that distinguishes the present 
analysis is that it applies to frequencies well below 
lu* where the scattering of transverse sound waves is 
weak, whereas I* was originally introduced to explain 
the plateau that exists in the density of states above lo* . 
Recent investigations indicate that the vibrational modes 
above lo* are diffusive and the corresponding mean free 
path is controlled by a distinct length scale id ~ At/) - * 
(see Sec. |VJ). 

IV. NUMERICAL TESTS 

In the previous section, we obtained the packing frac- 
tion dependence of the (inverse) attenuation coefficient in 
the Rayleigh scattering regime by introducing the char- 
acteristic length t* by hand. In order to lend some sup- 
port to this phenomenological approach, we show that 
the final result in Eq. (fT2]) is consistent with numerical 
evaluations of the Kubo formula for the energy diffusivity 
that do not rely on any ad hoc assumptions. 

A. Energy diffusivity 

In the weak scattering regime oj -c lo*, where 
Rayleigh's law holds, the energy diffusivity, d(uj), can 
be approximated as 

d(w) ~ vtts(u) (13) 

Inspection of Eq. (fl2l) and Eq. ((4]) allows us to write 

d(u) lu 4 ~ A# . (14) 

In Ref . [l| , the energy diffusivity of an unstressed pack- 
ing [34j], comprised of a 50/50 bidisperse mixture of 2000 
frictionless spheres with a diameter ratio of 1.4, was com- 
puted numerically using the Kubo formalism [15[ 

= I E(^«)~ 2 l4f 8{e>i ~ uj), (15) 
j 



where Sij denotes the energy flux matrix elements Sij. 
They read pj| 

g = (co ^ (f m -? n)ei (rn-a) fl$* e>;/3) . 

4 UJi OJj — 

7nn,ap 

(16) 

where r m indicate the particle positions, H™a the dy- 
namical matrix and ej(m; a) the corresponding eigenvec- 
tors [13]. 

The numerical evaluation of Eq. (|15|) carried out 
in Ref. [l[ provided clear of evidence for a sharp en- 
ergy transport crossover at a characteristic frequency 
lo* ~ A05 that separates the Rayleigh law regime from 
a diffusivity plateau. However, the packing fraction de- 
pendence of the pre-factor of Rayleigh law was not in- 
vestigated. We now proceed to re-analyze the numerical 
data obtained in Ref. [l[ to extract and test the critical 
exponent predicted by Eq. (fT4"|) for the diffusivity in the 
low lo regime . 



B. Results 

The scaling analysis of the previous section was implic- 
itly carried out in an infinitely large system for which vi- 
brational modes are continuously distributed in the spec- 
trum. However, in a cubic simulation box of size L the 
low lo acoustic waves occur close to the discrete frequen- 
cies allowed by the linear dispersion : 

Ui = ?ZpVp2 + (? 2 + r 2 ( 17 ) 
Li 

where {p, q, r} denote the Miller indexes and Vt the speed 
of sound. This allows the calculation of the mode degen- 
eracies predicted by classical elasticity For example 
the lowest lying mode (±1,0,0) should be 12- fold de- 
generate because there are two transverse directions for 
the six allowed wave vector of magnitude 2n/L. Simi- 
larly the second lowest lying mode would have a degen- 
eracy number equal to 24. In finite disordered samples 
these degeneracies are partially lifted even in the regime 
lo <C lo* where the vibrational modes are indeed quasi 
plane waves. 

The correct procedure to bin the data, is to average 
the diffusivity and frequency of each degenerate group 
of long wavelength phonons separately for each packing 
fraction considered. Figure [1] shows a log- log plot (based 
on the data published in Fig. 1 of Ref. [l| ) of d(uS) uj 4 
versus packing fraction A0 for the lowest lying group of 
12 modes at A</> = 0.3 , 0.1 , 0.05 , 0.02 and 0.01 (circles) 
and the second group of 24 modes at A0 = 0.1 , 0.05 
and 0.02 (squares). The continuous line was drawn with 
slope 11/4, as predicted in Eq. (TT4]) . This scaling appears 
to be consistent with the numerical results presented in 
Fig. [TJ thus providing indirect evidence for the length 
scale I* . Larger system sizes are needed to further test 
the validity of our findings over a greater range of A<fi. 
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FIG. 1: Log-Log plot of d(u) uj 4 versus packing fraction A(j> for 
the lowest lying group of 12-fold degenerate modes at A(f> = 
0.3 , 0.1 , 0.05 , 0.02 and 0.01 (circles) and for the second 
lowest lying group of 24-fold degenerate modes at A<f> = 0.1 , 
0.05 , 0.02 (squares). The slope of the continuous line is 11/4. 
The data is reproduced from Fig. 1 of Ref. [1J. 



The main requirement is that the packing fraction is large 
enough to ensure that the plane wave frequencies are well 
below uj* ~ (4> — ^c) 5 , but still sufficiently close to cf) c to 
be in the critical regime. 



ONSET OF THE STRONG SCATTERING 
REGIME 



Above uj*, the scaling analysis of Eq. (|T4|) does not 
apply because the diffusivity plateaus to a minimal value 
do independent of packing fraction [l|, 0] . This plateau 
starts at uj* and ceases to exist only at the very end of 
the spectrum where localization sets in and the diffusivity 
becomes vanishingly small. 

We will now provide a brief summary of a scaling anal- 
ysis detailed in Ref. [l| for the packing fraction depen- 
dence of uj* that marks the onset of the plateaus in the 
diffusivity and density of states. Since do is weakly de- 
pendent on Acj) the mean free path at the onset of the 
diffusivity plateau is given by |l| 



do 



A<j> 



1/4 



(18) 



where we used Eq. (fT3| and Eq. (j4j. The diverging 
length scale £d that controls the onset of the diffusive 
modes above uj* is distinct and larger than £*, for small 
Acf). 

Although id can be arbitrarily large compared to the 
particle size, it is comparable to the corresponding wave- 
length A* at uj* . As a result, energy transport is severely 



hindered above uj* . The validity of this heuristic crite- 
rion for the onset of strong scattering, l d ~ A*, has been 
corroborated by numerical tests in Ref. 0, 0| . This al- 
lows one to readily convert the wavelength (of order ij) 
to the corresponding frequency uj* by using the phonon 
dispersion relation for transverse phonons. The result 
reads 



(19) 



A natural question that arises in this context is how 
the low frequency regime of diverging diffusivity " meets" 
the plateau in d(uj) expected above uj* . The plot of the 
diffusivity obtained from the numerical evaluation of the 
Kubo formula in Ref. [l|, Hf lacks the spectral resolution 
necessary to probe in detail the crossover region. Recent 
mean field calculations carried out on lattice models, that 
mimic several of the distinctive properties of soft spheres 
packings, predict that the pre-factor of Rayleigh's law 
scales as A(f> 7 ^ 4 for uj <C uj* [Uj]. This result differs from 
the Acf) 11 / 4 scaling obtained in Eq. (fl~4|) . This discrep- 
ancy may originate from the fact that the mean field 
model of Ref. 32| does not capture the length scale £* 



which is the crucial ingredient that enters our scaling ar- 
gument leading to Eq. (fl~4|) . 

In the mean field calculations, the Rayleigh scattering 
regime extends almost all the way to frequencies of order 
uj* so that cLmf{u*), the mean field diffusivity at uj* , is 
given by 



IMF 



{uj*) 



A<t> 



7/4 



1 



UJ" 



A(f>- 



(20) 



This implies that, as the frequency is lowered through 
uj* , the diffusivity exhibits an upward jump of magnitude 
A<j)~4 which becomes arbitrarily large as <f> approaches 
4> c . By contrast, dEx(oJ*), the value of the diffusivity 
extrapolated to uj* according to Eq. (fLrj) reads 



(uj*) 



A<j) 



11/4 



Acj)^ 



(21) 



Note that there is no reason why the Rayleigh regime 
ought to extend all the way to w*. As this critical fre- 
quency is approached, the scattering becomes stronger 
and the assumptions under which Eq. (|14j) was derived 
no longer hold. Nonetheless, if such an extrapolation is 
made, the intercept of the low uj Rayleigh branch of the 
diffusivity will fall below the plateau value do, eventually 
becoming vanishingly small as <p approaches <fi c . This sce- 
nario would imply a non-monotonic uj dependence of the 
diffusivity analogous to the behavior of the participation 
ratio 331. 



VI. 



CONCLUSIONS 



In this article, we have shown that the elastic atten- 
uation coefficient a(uj) of long-wavelength shear sound 
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waves exhibits critical scaling near point J according to 
the power law a(u) ~ w 4 (</> — </> c )~ 5 / 2 and diverges as <f> 
approaches <j) c . A simple scaling argument has been pro- 
posed to obtain the packing fraction dependence of the 
attenuation coefficient which relies on two crucial physi- 
cal ingredients: the scaling of the shear modulus G that 
vanishes at point J and the existence of a characteris- 
tic length scale I* that characterizes the spatial fluctua- 
tions in the structure and diverges at = <f> c . The crit- 
ical exponent suggested in this work is consistent with 
numerical evaluations of the Kubo formula for the en- 
ergy diffusivity at low u>. The favorable comparison be- 
tween our scaling argument and numerics supports the 
simple view that a jammed solid can be regarded as a 
mosaic composed of domains of characteristic size given 
by I* ~ (cj) — (f) c )~2. We hope that the effects discussed in 



this work will stimulate further experimental investiga- 
tions of the dependence of sound attenuation in emulsions 
and granular media under loading. 



Acknowledgments 

I thank N. Xu, M. Wyart, A. J. Liu and S. R. Nagel 
for helpful discussions. I am grateful to B. Tighe, M. van 
Hecke, Y. Shokef and W. T. M. Irvine for their comments. 
This work was initiated during my stay at the Weizmann 
Institute under the auspices of the Feinberg Foundation 
Visiting Faculty Program Fellowship which I gratefully 
acknowledge. 



[1] V. Vitelli, N. Xu, M. Wyart, A. J. Liu and S. R. Nagel, 

Phys. Rev. E 81, 021301 (2010). 
[2] N. Xu, V. Vitelli, M. Wyart, A. J. Liu, and S. R. Nagel, 

Phys. Rev. Lett. 102, 038001 (2009). 
[3] Physics of Amorphous Materials., S. R. Elliott, Long- 
man, Harlow (1990), pp. 222-243. 
[4] F. Leonforte, R. Boissiere, A. Tanguy, J. P. Wittmer and 

J. -L. Barrat, Phys. Rev. B 72, 224206 (2005). 
[5] M. Wyart, Annales de Phys. 30 (3), 1 (2005). 
[6] V.L. Gurevich, D.A. Parshin and HR. Schober, Phys. 

Rev. B 71, 014209 (2005) and references therein. 
[7] Amorphous Solids. Low Temperature Properties, edited 

by W. A. Phillips, Springer- Verlag, Berlin (1981). 
[8] X. Jia, C. Caroli, B. Velicky, Phys. Rev. Lett 82, 9, 

(1999). 

[9] C.H. Liu and S. R. Nagel Phys. Rev. Lett. 68, 2301, 

(1992) . 

[10] X. Jia, Phys. Rev. Lett. 93, 154303 (2004). 
[11] T. Brunet, X. Jia and P. Mills, Phys. Rev. Lett. 101, 
138001 (2008). 

[12] T. Brunet, X. Jia and P. A. Johnson, Geophys. Res. Lett. 

35, L19308 (2008). 
[13] V. Lubchenko and P. G. Wolynes, Proc. Natl Acad. Set. 

USA 100, 1515 (2003). 
[14] P. Sheng and M. Zhou, Science 253, 539 (1991). 
[15] P. B. Allen and J. L. Feldman, Phys. Rev. B 48, 12581 

(1993) ; J. L. Feldman, P. B. Allen, and S. R. Bickam, 
Phys. Rev. B 59, 3551 (1999). 

[16] D. G. Cahill and R. O. Pohl, Phys. Rev. B 35, 4067, 
(1987). 

[17] X. Liu and H. v. Lhneysen, Europhys. Lett. 33, (8), 617 
(1996). 

[18] W. Schumacher, G. Ruocco and T. Scopigno, Phys. Rev. 

Lett. 98, 025501 (2007). 
[19] D. A. Parshin and C. Laermans, Phys. Rev. B63, 132203 

(2001). 



[20] M. Wyart, S. R. Nagel, and T. A. Witten, Europhys. 

Lett. 72, 486 (2005). 
[21] M. Wyart, L. E. Silbert, S. R. Nagel, and T. A. Witten, 

Phys. Rev. E72, 051306 (2005). 
[22] J. H. Page, W. K. Hildebrand, J. Beck, R. Holmes, and 

J. Bobowski Phys. Stat. Sol 1, 11, 2925 (2005). 
[23] J. E. Graebner, B. Golding, and L. C. Allen, Phys. Rev. 

B, 34, 5696 (1986). 
[24] C. S. O'Hern, L. E. Silbert, A. J. Liu, and S. R. Nagel, 

Phys. Rev. E 68, 01136 (2003). 
[25] L. E. Silbert, A. J. Liu, and S. R. Nagel, Phys. Rev. Lett. 

95, 098301 (2005). 
[26] B. J. West and M. F. Shlesinger, J. of Stat. Phys., 36, 

5, (1984). 

[27] S. John, H. Sompolinsky, and M. J. Stephen, Phys. Rev. 
B 27, 5592 (1983). 

[28] D. J. Durian, Phys. Rev. Lett. 75, 4780 (1995). 

[29] S. Alexander, Phys. Rep. 296, 65 (1998). 

[30] N. W. Ashcroft and D. N. Mermin, Solid State Physics, 
Brooks Cole, (1976). 

[31] W. G. Ellenbroek, E. Somfai, M. van Hecke, W. van 
Saarloos, Phys. Rev. Lett. 97, 258001 (2006). 

[32] M. Wyart, Europhys. Lett. 89, 64001 (2010). 

[33] N. Xu, V. Vi telli, A. J. Liu, and S. R. Nagel, 
larXiv:0909.3701l 

[34] The unstressed packings are obtained by replacing the in- 
teraction potential, V(r;j), between each pair of overlap- 
ping particles with an unstretched spring with the same 
stiffness, V"(r'?), where r\^ is the equilibrium distance 
between particles i and j. This corresponds to dropping 
terms depending on the first spatial derivative of the po- 
tential, V' , in the dynamical matrix which increases the 
mode energy and frequency, because V' is negative for 
repulsive interactions. 



